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ABSTRACT 

We present an analysis of Chandra observations of the central regions of the cooling flow cluster Abell 
2029. We find a number of X-ray filaments in the central 40 kpc, some of which appear to be associated 
with the currently active central radio galaxy. The outer southern lobe of the steep-spectrum radio 
source appears to be surrounded by a region of cool gas and is at least partially surrounded by a 
bright X-ray rim similar to that seen around radio sources in the cores of other cooling fiow clusters. 
Spectroscopic fits show that the overall cluster emission is best fitted by either a two temperature 
gas {kThigh = 7.47 keV, kTiow = 0.11 keV), or a cooling flow model with gas cooling over the same 
temperature range. This large range of temperatures (over a factor of 50) is relatively unique to 
Abell 2029 and may suggest that this system is a very young cooling flow where the gas has only 
recently started cooling to low temperatures. The cooling flow model gives a mass deposition rate of 
M = 56^21 -^0 yr~^- In general, the cluster emission is elongated along a position angle of 22° with 
an ellipticity of 0.26. The distribution of the X-ray emission in the central region of the cluster is 
asymmetric, however, with excess emission to the north-east and south-east compared to the south- 
west and north-west, respectively. Fitting and subtracting a smooth elliptical model from the X-ray 
data reveals a dipolar spiral excess extending in a clockwise direction from the cluster core to radii of 
~ 150 kpc. We estimate a total mass of Mspirai ~ 6 x 10^^ Mq in the spiral excess. The most likely 
origins of the excess are either stripping of gas from a galaxy group or bare dark matter potential 
which has fallen into the cluster, or sloshing motions in the cluster core induced by a past merger. 
Subject headings: cooling flows — galaxies: clusters: general — galaxies: clusters: individual 

(Abell 2029) — intergalactic medium — radio continuum: galaxies — X-rays: 

galaxies: clusters 



Ph' 1. INTRODUCTION 

Q ■ An intracluster medium (ICM) of thermal X-ray gas 
^ fills galaxy clusters and traces the structure within the 
c/) . cluster's gravitational potential. In a relaxed system, 
■ the thermal gas in outer regions of the cluster has a rel- 
^ ' atively smooth and symmetric surface brightness distri- 
• ^ • bution, while the inner regions often reveal peaked X-ray 
' emission which is generally interpreted as a cooling flow 
5h (Fabian 1994). The thermal gas in the central regions 
. . 1 of these cooling flow clusters has a radiative cooling time 
which is significantly shorter than the Hubble time, lead- 
ing to mass deposition rates of more than 100 Mq yr~^ 
being determ ined from imag ing data from Einstein and 
ROSAT fe.g. lWhite e'ralll9 97: Allen 2000). Analysis of 
ASCA SIS and CIS data su ggested muc h lower inferred 
cooling rates f^seelMakishim aet al.ll200'lL and references 
therein), and more recently, high-resolution spatially re- 
solved spectral fitting of Chandra and XMM data have 
revealed cooling rates that are a factor of a few to ten 
times l ower than those o btained from the imaging analy- 
sis (^see lMcNamaral2002l and references therein). Further 
evidence of lower cooling fiow rates comes from obser- 
vations using XMM's Refiection-Grating Spectrometer 
which finds no evidence of cool gas in cluster cores below 
approximately one-third of the ambient cluster temper- 
ature (e.g., Peterson et al. 2003). 

The cores of cooling flow clusters typically host giant 
cD galaxies which often contain an excess blue stellar 
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component, indicating t hat there is recent star formatio n 
in these central regions ljMcNa,maTa fc 0'Connel]HT98l . 
The cores of these clusters also often contain powerful, 
steep spectrum radio galaxies. Chandra observations 
of these systems reveal complex surface brightness dis- 
tributions with X-ray depressions surrounded by bright 
rims (e.g., Perseus, Fabian et al. 2000; Hydra A, Mc- 
Namara et al. 2000; AbeU 2052, Blanton et al. 2001). 
The X-ray depressions are frequently filled with radio 
synchrotron plasma from the lobes of t he centra l radio 
source. The theoretical models of iHeinz et all 1119981) 
predicted that the expanding radio source would drive 
shocks into the ICM and heat the surrounding gas. In 
contrast to this model. X-ray observations reveal that the 
rims surroundi ng the d epression are composed of cool 
gas (e.g., Sc hmidt et al . 2002: NTflsen et al.ll200l . al- 
though more recent models suggest that the rims may 
be for med by weak shocks d riven by the expanding radio 
lobes ijRevnolds et alJl200H) . In addition, X -ray observa- 
tions of the Perseus cluster by^bian et al. (2003) may 
provide evidence of weak shocks or sound waves driven 
into the ICM by the central radio source. Estimates of 
the total energy output from the central radio source 
show that in at least some cases it is sufficient to off- 
set the effects of the cooling flow (e.g., Virgo A, Owen, 
Eilek, & Kassim 2000; Hydra A, David et al. 2001; AbeU 
2052, Blanton, Sarazin, & McNamara 2003, Birzan et al. 
2004). 

In this paper we present results of a radio and X- 
^q^gamafc^s of the cooling flow cluster Abell 2029. 
Abell 2029 is a nearby, z=0.0767, Bautz-Morgan type 
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I, richness class 4.4 l|Dressleilll978j) cluster of galaxies. 
Optical observations show that diffuse light from the 
central cD galaxy (IC 1101 ) extends over more than 
600 kpc (Uson, Bo ughn. fc K uhn 1993i), making it one 
of the largest known galaxies. The central region of 
Abell 2029 is host to the steep spectrum radio source 
PKS 1508+059. The inner regions of the radio galaxy 
show a compact flat-spectrum core and two oppositely di- 
rected jets which undergo roug hly 90° bends at dis tances 
of ~ 10 - 15" from the core ijTavlor et all 119941. Be- 
yond the bends. iTay lor et aP find that the radio emission 
extends into two steep-spectrum radio tails. The over- 
all source morphology is a C-shaped wide-angle-tailed 
(WAT) s tructure that is typical of cluste rs undergoing 
mergers l|Roettiger. Burns, fc Lo kcn 1996), although the 
source size of 35 kpc is smaller t han average for WATs 
llO'Donoghue. Eilek. fc (>^IT99(lD . 

The intracluster medium (ICM) in Abell 2029 has 
been extensively studied in X-rays. On large scales, 
the X-ray emission from the cluster appears to be very 
relaxed (Buotc fc Canizarcs 1996; _Sarazin ct al. 199S 
iMolendi fc De Gran di 1999), consistent with the large 
inferred cooling flow rates of M = 2 00 - SOOM^yr"^ 
llSarazin et all 1992: Edge e^]^9^ iPeres et all 199S ; 
Sarazin et all Il99a although see also lWhiteri200C : 
Lewis et all I2002D Optical observations are somewhat 



inconsistent with the relaxed cooling- flow picture as they 
do not reveal any of the typical star-formation indi- 
cators s uch as O II or blue stellar c olors in the cen- 
tral cD fMcN amara fc O'ConneU * 1989). In this respect, 
Abell 2029 is nearly unique among clusters with very 
large inferred cooling rates. 

Higher res olution ROSAT HRI observations by 
iSarazin et al.l |l992) of the core reveal a number of X-ray 
filaments in the central 30" (43 kpc). iTavlor et all l|1994l) 
find that the radio emission appears to be anti-correlated 
with the filaments, suggesting that the radio structure 
traces the low X-ray pressure regions. The presence of 
the X-ray filaments is still controversial, ho wever, as sub- 
sequent analysis of the ROSAT HRI data bv I White et al.l 
((i994) found no significant evidence of filaments. 

Throughout t his paper we adopt WMAP cosmologi- 
cal parameters l|Bennett et al.ll200^ Hq = 71 km s~^ 
Mpc-\ f^A=0.73, and n„^=0.27. At the redshift of Abell 
2029, this corresponds to a scale of 1.44 kpc/arcsec. The 
uncertainties quoted in this paper are 90% confidence 
intervals. 

2. OBSERVATIONS AND DATA REDUCTIONS 

Abell 2029 was observed with Chandra on 2000 April 
12 for a total of 19.8 ksec. The observations were 
taken in Faint (F) mode with the cluster center on the 
S3 chip which was operating at —120 C. We extracted 
the archival observations (OBSID 891) and reprocessed 
the data using CIAO v2.3, and CALDB v2.18 (except 
where noted otherwise). Only events with ASCA grades 
0,2,3,4, and 6 were included in the analysis. The back- 
illuminated SI CCD chip was used to examine the back- 
ground during the observations as the cluster emission 
fills the entire S3 chip. The data were processed using 
the lc_clean script of M. Markevitch to provide consis- 
tent comparison with the background files. The data 
were free of large flares, and we have removed only 128 s 
of exposure. The period C ACIS blank-sky background 



data files included in the Chandra Calibration Database^ 
were used in the data analysis. 

3. X-RAY IMAGE 

The raw Chandra image of the entire S3 chip (roughly 
860 X 860 kpc) is shown in Figure ^ This image covers 
the 0.3 — 10 keV energy band and has not been corrected 
for exposure or background but has been binned by two 
pixels. The broad cluster core is visible at the center of 
the image, and is surrounded by more extended diffuse 
X-ray emission. The cluster emission is extended in the 
north-east to south-west direction, and is not distributed 
symmetrically about the cluster core. In particular, there 
is a strong decrease in surface brightness ~ 20 kpc to the 
west of the cluster core and a much more gradual de- 
crease in the other directions. More generally, there ap- 
pears to be excess emission to the north-east and south- 
east compared to the south-west and north-west, respec- 
tively. This lack of mirror-symmetry has been previously 
reported from the ROSAT HRI data l|Sarazin et al.ll992l) 
and is discussed in more detail in §E1 
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Fig. 1.— Raw Chandra image of the full ACIS S3 chip (~ 860 X 
860 kpc) region of Abell 2029 in the 0.3 - 10.0 keV energy band. 
The image has not been corrected for background or exposure but 
has been binned by a factor of two. The X-ray emission is elongated 
in the north-east to south-west direction and shows excess emission 
to the north-east and south-west of the core compared to the south- 
west and north-east. The linear negative east-west feature 1'.5 
below t he cluster center i s a shadow cast by a foreground spiral 
galaxy IClarke et aljl2(X)l) . 

We have identified individual X-ray sources in the 
central regions of Abell 2029 using the wavdetect al- 
gorithm in CIAO (Freeman. Kas hvap. Rosner. fc LambI 
[2002). The significance threshold was set at 10~^ which 
corresponds to roughly 4.7cr, or less than 1 false detection 
in the S3 field of view. Using the background for each 
source, we excluded two sources near the detector edges 
that had a signal-to- noise ratio of less than 3. The final 
list from the detection algorithm consists of 17 sources. 
These sources are shown overlaid on the DSS II red im- 

^ http://cxc.harvard.edu/caldb/ 
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age in Figure |2] In the central region of the cluster the 
algorithm finds both a peak around the cluster-center 
galaxy IC 1101, and one other peak which appears to 
be associated with structure in the X-ray filament which 
traces the radio tail (see § El . We have compa red the 
sourc e positions to the USNO A2.0 catalog (Mon etl;t al.l 
11998(1 and find 7 matches (including IC 1101) for coun- 
terparts within 2'.'8. Four of the optical counterparts fall 
less than 1" from the position of the X-ray peak found 
by wavdetect. Using the T wo Micron All Sky Survey 
catalog f2MASS: IC^tri et al, 2001) we find three sources 
which fall within 1" of the X-ray positions. Compar- 
ing the offsets of these sources, we find that there is no 
evidence of a position shift in the X-ray image, thus we 
have not applied additional astrometry corrections to the 
X-ray data. 



clearly asymmetric. 
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Fig. 2. — DSS II red image of the region covered by the Chandra 
S3 chip. The large object near the center of the image is the cluster- 
center cD galaxy IC 1101. The linear optical source 1.'5 south of 
the cD core is the edge-on spiral galaxy which produces photo- 
electric absorption of the X-ray emission (see [Clarke ct al. 200^. 
The X-ray sources identified by wavdetect with (without) 2MASS 
identifications are indicated in blue (red). 

Figure O shows an adaptively smoothed 0.3 — 10 keV 
Chandra image of the central ~ 1.'5 x lf5 region of 
Abell 2029. The adaptive smoothing was done using 
csmooth within CIAO and used a minimum S/N of 3 
per smoothing beam. The image has been corrected for 
background and exposure. This smoothed image shows 
more clearly the structure in the central regions of Abell 
2029. The broad cluster core displays an hourglass shape 
and there are a number of X-ray enhancements visible as 
seen in the raw data. On larger scales, the emission is 
extended in the north-east to south-west direction. The 
X-ray elongation is similar to that of the diffuse optical 
halo emission surroundin g the c entral cD. i?-Band ob- 
servations by lUson et al] 1(19911) trace the diffuse light 
to over 600 kpc from the core along a position angle of 
21° east of north. The ellipticity and position angle of 
the X-ray emission are discussed in more detail in § El 
However, unlike the optical emission, the X-ray image is 
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Fig. 3. — Adaptively smoothed X-ray image of the central ~ 
I'.b X I'.b region of Abell 2029 in the 0.3 - 10.0 keV energy band. 
The image has been corrected for the effects of background and 
exposure. The image shows the broad cluster core surrounded by 
a number of X-ray filaments which extend to the south-east, south- 
west, north-west, and north-east. 

Previous ROSAT HRI images of the inner 30" of Abell 
2029 showed a number of X-ray filaments running to the 
north -east, east, south, and south-west ijSarazin et alJ 
I1992() . Although the Chandra data do not precisely re- 
produce the locations of the filaments. Figure 13 shows 
clear X-ray enhancements to the north-east, east, south- 
west, and north-west of the core. If we just consider 
the south filament tracing the radio lobe, which we con- 
centrate on in this paper, we find that the counts in a 
rectangle surrounding the filament are more than 9(t in 
excess of the average counts in an annulus centered on 
the core and covering the range of radii of the filament. 

4. CLUSTER PROPERTIES 
4.1. Integrated X-ray Spectrum of the Cluster Center 

Our goal in this section is to obtain the best spectral 
fit to the integrated spectrum of the cooling fiow region 
of Abell 2029. We begin our analysis with the simplest 
model (single temperature MEKAL) before proceeding to 
the more complex models. These spectral fits differ from 
those of Lewis ct al. (2002) who investigated the radial 
distribution of spectral parameters in concentric annuli, 
rather than the integrated spectral properties. 

The spectral properties of the central region of Abell 
2029 were studied using the XSPEC vll.2 software pack- 
age (pVrnaud 1996). The 0.7 — 8.0 keV spectrum was ex- 
tracted for a circular region of radius 116" ( 167 kpc), cor- 
respon ding to the cooling radius found bv lSarazin et all 
l(1992j) . centered on the cluster core. The extracted re- 
gion excluded point sources discussed in §|2|but retained 
the core as no obvious power-law AGN contribution is 
seen (see § 15. 1|) . The CIAO tool ac is spec was used 
to extract the spectrum and produce the weighted re- 
sponse files. The spectrum contained a total of 163,000 
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background-subtracted counts. We used the more recent 
CALDB v2.26 under CIAO v3.0.2 for the spectral ex- 
traction as it incorporates the best model of the effects 
of the ACIS QE degradation in the response file. The 
spectrum was binned in energy using the grppha task so 
that each bin contains at least 25 counts. The blank- 
sky background files discussed in § |21 were extracted in 
the same manner over the circular region, excluding the 
regions around sources detected in the target field. The 
spectral fits showed significant residuals in the ^ f .4— 2.2 
keV energy range that are likely due mainly to the ACIS 
calibration problems around the Si-K and mirror Ir fea- 
tures, but may also include a contribution from ICM Si 
and S lines. To exclude the majority of the residuals we 
would need to exclude this entire energy range. Our fits 
showed that the addition of a systematic error term of 
3% to the data allowed us to retain the majority of the 
spectral region in the fits. The spectral fits (including the 
systematic error) were done using both the full 0.7 — 8.0 
keV band (summarized below) and also with the 1.8 — 2.1 
keV region around the mirror Ir edge excluded. We note 
that fits excluding the 1.4—2.2 keV regions without the 
addition of a systematic error term are consistent with 
those presented here. The results of the spectral fittin g 
from both energy ranges are summarized in Table fSTI 
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Fig. 4. — Spectrum of the central 116" cooling core region. The 
spectrum is fitted by two MEKAL models with a systematic error term 
of 3% added to the data. The absorption was fixed to Galactic and 
the metallicity of the second model was fixed to the (fitted) value 
of the first model. 

Our initial spectral fitting of the central region of Abell 
2029 used the single temperature MEKAL plasma em ission 
code (|Kaastra fc Mewelll993t iLiedahl et aLlllQQRD to fit 
the spectrum in the energy range of 0.7 — 8.0 keV. Fix- 
ing the absorption to the Galac tic value of 3 . 14 x 10^° 
cm~^ as determined from Dickev fc LockmanI ljl99ClD re- 
sults in a reduced of 1.23, and yields a temperature 
of kT ^ 7.27t:i^ keV and an abundance of Z = 0.47t;^^ 
times solar (assuming metallicity of And ers fc Grevessd 
[I989i) . Given the relatively poor quality of this fit, we 
have allowed the absorption to be a free parameter and 
find that the resulting model is somewhat improved (re- 
duced = 1.18), but this fit requires the absorption 
column to fall below Galactic which may be unlikely. It 
is more likely that this fit indicates that there is an addi- 
tional soft X-ray emission component (eg. a second low 



temperature thermal component) which is not modeled 
by a single temperature, and/or the QE degradation cor- 
rection may not be exact. 

The spectral fits are improved (at > 99.9 % level ac- 
cording to the F-test) by adding a second MEKAL com- 
ponent to the model. In this model, the abundances of 
the two MEKAL components are tied together as there is 
not enough information in the spectrum to accurately 
fit them individually. For the two temperature fit with 
absorption fixed to Galactic, we find a reduced of 
1.18. The best fitting high and low temperatures are 
fcThigh = 7A7tf5 keV and fcTiow - 0.1ltEJ°^ keV, 
and the abundance is Z = 0.49^'^ Qg times solar. The 
emission measures (MEKAL normalization parameter K) 
of the hot and cold components are within a factor of 
about two of each other [Khigh = 5.1 x 10~^ cm~^ 
and Kiow — 2.1 x 10^^ cm ~^). Our metallicit y results 
are in good agreemen t with i Lcwis et al.' ("2002) and the 
BeyyoSAX res ults ofl Molendi fc De Grandi (1999 ) and 
llrwin fc Bregm an (2001). This best-fit two-temperature 
model is shown in Figure 01 Allowing the absorption 
to be free results in a model with the same reduced x^j 
temperatures and abundance to within the 90% confi- 
dence errors, and only slightly lower than Galactic ab- 
sorption of Nh = 2.22;;75 x 10^° cm'^. We find that 
the temperatures of the hot and cool components of the 
two MEKAL fits differ by more than a factor of 50. One 
possible explanation for the soft emission component in 
Abell 2029 is emission from filaments of W arm-Hot Inter- 
galact ic Medium surrounding the cluster. iKaastra et alJ 
1)2003(1 have recently studied a sample of 14 clusters us- 
ing XMM-Newton data and find evidence for a spatially 
extended soft component in the spectrum of five clusters. 
They find that the emission measure of the soft excess 
is limited to roughly 10% of that of the hot gas compo- 
nent in the central regions of the clusters but increases 
significantly to the outer regions of the clusters, where 
the soft excess becomes dominant in three of the clus- 
ters. The small field of view of the ACIS-S3 ccd does 
not allow us to examine the soft component beyond the 
central regions of the cluster. We note, however, that the 
emission measure of the soft component in Abell 2029 is 
a much larger percentage (roughly 50%) of that of the 
hot component. This much larger central concentration 
in Abell 2029 suggests that the emission is not primarily 
due to intercluster filaments. 

We have considered the effects of uncertainties in the 
Chandra calibration and the Galactic absorbing column 
on the spectral fits since both could have an impact on 
the fitted soft X-ray component. The observations were 
taken early in the Chandra mission, thus the effects of 
the QE degradation uncertainty are small. To further 
limit any uncertainties from calibration problems at low 
energies we have restricted our spectral analysis to ener- 
gies above 0.7 keV. In the case of M87, high resolution H I 
observations showed Galactic col umns were r oughly 28% 
lower toward the cluster center l)Lieu et alJ ll996') than 
those determined by the lower spatial resolution observa- 
tions of Stark et al. ( 1993). We have ex amined the range 
of Nh measurements from iDickev fc LockmanI l|199r)D in 
the vicinity of Abell 2029 and find that they vary by at 
most 12% from the Galactic column used above. Fitting 
the spectra with a 12% (or even 30%) lower Galactic 



The Complex Cooling Core of Abell 2029 



5 



column does not change the spectral results within the 
90% confidence interval, ft is therefore unlikely that the 
soft emission seen in Abell 2029 is a result of either a 
calibration error or a variation in the Galactic column. 

Given the indications that the cluster has cool gas 
in the core, and the past high cooling flow rates of 
M > 200 Mq yr~i determined for AbeU 2029, we have 
also investigated a cooling flow model for the core. The 
spectral fits used the MKCFLDW model in conjunction with 
a MEKAL model such that the high temperature compo- 
nent (fcThigh) of the coohng flow model is tied to the 
temperature of the MEKAL model. The MEKAL model 
represents the hot outer cluster gas along the line of 
sight. We assume that the cooling flow material cooled 
from this ambient cluster gas, and thus we fix the abun- 
dance in MKCFLOW to that of the MEKAL model for the 
ambient cluster gas. When the absorption was fixed to 
Galactic, the MEKAL + MKCFLOW model had best fit high 
and low temperatures of /cThigh ~ '^■^'^-0 23 ^^'^ 
kTio^ = OMtoil keV, an abundance of Z = O.bOtom 
times solar, and a mass deposition rate of M — 56^21 -^^0 
yr~^ (reduced — 1.19). Note that the value of /cTiow is 
consistent with zero and the best-fit value is well below 
the lowest photon energies in the spectrum; thus, this fit 
represents gas which is cooling down to very low temper- 
atures. The data are not sufficient to constrain the MEKAL 
+ MKCFLOW model with the absorbing column, metallic- 
ity, kThigh, kTiow and cooling flow rates free. We have 
therefore fixed the lower temperature component in the 
cooling flow model to 0.08 keV based on the results from 
the MEKAL fits. The subsequent MEKAL + MKCFLOW fits 
have reduced = 1-18, a mass deposition rate of zero 
and an (unlikely) absorbing column ~ 4 times lower than 
Galactic. The upper and lower ranges of the tempera- 
ture distribution spanned by the cooling fiow model are 
consistent with the hot and cool components of the two 
temperature fits. This cooling flow temperature range 
of kThinh / kTinv, > 50 is far g reater than seen for other 
clusters (jPeterson et al.l 120031) . and is discussed further 
in §13 

We have also considered an annulus centered on the 
cluster core running from the cooling radius out to the 
edge of the S3 chip (i.e. 116" to 174", or 167 to 251 
kpc). We have extracted the spectrum for this region 
and examined both single and two temperature MEKAL 
models as well as the MEKAL + MKCFLOW model. The 
results of these fits are shown in Table IA2I As with 
the cooling radius fits, we find that the two tempera- 
ture model and the MEKAL + MKCFLOW provide better fits 
than the single temperature model (i.e. x^/d.o.f.MBKAL 
= 436/343, xVd-o.f ■MEKAL+MEKAL = 399/341, and 
x'^/d.oI.MEKAL+MKCFLOW=^0Q/^4:l). The high tem- 
perature component from the best fit models in this 
outer annulus is somewhat higher than that for the cool- 
ing core, consistent w i th the radial temperature increase 
found bv 'Lewi s et al.l l)2002() . The low temperature com- 
ponent remains around 0.1 keV indicating that the cool 
gas continues to radii of at least 250 kpc from the core 
of AbeU 2029. 

The multi-temperature gas and moderate cooling rate 
we have found for the spectral fits to the Chandra data 
of AbeU 2029 are in a possible conflict with the analysis 
of iLewis et al.l (l|2002l) who found no evidence of multi- 



ph ase gas in th e cluster. The spectral fits undertaken 
bv lLewis et all covered the energy range of 0.3 — 8.0 keV 
but did not include the (then unknown) correction for the 
QE degradation at soft energies. Without the QE cor- 
rection, the spectral models would over-predict the flux 
at low energies thus resulting in weaker constraints on a 
cool spectral component. The discr epancy between our 
multi-temperature fits and those of ILewis' "eTaLl ll200l 
are likely due in part to the improved calibration at low 
energies provided by CALDB v2.26 as well as our re- 
striction to energies above 0.7 keV. We no te also that 
the ni ultiphase and coohng flow models of ILewis et alJ 
l)2002l) were only applied to each of their inner three an- 
nuli which cover radii < 30", while our models were ap- 
plied to a single region out to the cooling radius of 116". 

Based on the models we examined, we find that the 
best fit is produced by the two-temperature MEKAL model 
where the cluste r ga s has kThigh = 7.47 keV and kTiow = 
0.11 keV (TablelMj. This fit is only slightly better than 
the MEKAL + MKCFLOW model with absorption fixed to 
Galactic where we find gas cooling over a range of tem- 
peratures between 7.94 keV and 0.08 keV, and a mass 
deposition rate of M = 56^^^ Mq yr~^. 
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Fig. 5.— The central 50" X 50" (72 x 72 kpc) of the adaptively 
smoothed Chandra image shown in Figure jSJ The c ontours show 
the 14 90 MHz radio emission of PKS 1508-1-059 from lTavlor et alJ 
11994) . There are several filaments connected to the currently ac- 
tive radio source: the southern radio lobe is traced by an X-ray 
filament to the east, and the northern jet is also traced by a fila- 
ment. There are also several filaments to the north-cast which are 
apparently unrelated to the current radio source. 

5. CLUSTER CENTER X-RAY AND RADIO INTERACTION 

The central region of Abell 2029 is host to the 
C— shaped radio source PKS 1508-1-059. T he central ra- 
dio source has been extensively studied bv lTavlor et alJ 
{^94). In Figure we show the 1490 MHz radio con- 
tours from Taylor et al. overlaid on the central 50" region 
of the adaptively smoothed Chandra image. The radio 
source is composed of two inner jets which extend out to 
distances of ^ 10 — 15" and two outer tails with sizes of 
^ 20—25" 30—35 kpc). Both radio tails are connected 
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to the inner radio jets by faint synchrotron bridges. The 
southern radio tail is bounded to the north by an X-ray 
filament which traces the full length of the tail. A simi- 
lar correspo ndence betwee n the radio and X-ray emission 
was seen bv lTavlor et al .' when they compared the 1490 
MHz radio data to the ROSAT HRI image. 

The thermal pressure in the X-ray filament was de- 
termined by extracting a rectangular region set on the 
filament and fitting the data within XSPEC. The effect 
of the overlying ICM was modeled by using a region just 
exterior to the filament as a local background during the 
fits. The spectra were binned to 25 counts/bin and were 
fitted with a single temperature MEKAL model. The ab- 
sorption was fixed to Galactic, and the abundance was set 
io Z = 0.45 times Solar. The best-fit model {x^ /d.o.f. = 
38/40) gave a temperature of kT = A.2At]:%, keV, and 
a normalization oi K = 1.62 ± 0.20 x lO^'' cm^^, where 
the normalization is given by 

Here, Da is the angular diameter distance, z is the red- 
shift. He and uh are the electron and hydrogen densities 
in cm~^, respectively. We assume that the emission is 
from a prolate cylinder and we find an electron density 
of rig = 0.10 cm~'^ and a pressure of Pth = 1.6 x 10~^ 
dyn cm~^ . Using the same XSPEC model and local back- 
ground we fitted the emission on either side of the fila- 
ment and find that the filament is in pressure equilib- 
rium with the ambient thermal gas. The thermal pres- 
sure is a factor of 50 larger than the minimum energy 
synchrotron press ure of Pms = 2.7 x 10~^^ dyn cm~^ 
as determined bv iTavlor et al.1 ljl994|) and converted to 
our cosmology. Using the minimum energy parameters 
from^avlor et al. (1994), the ratio of the energy in the 
relativistic protons compared to that in electrons (fc), 
to the volume filling factor of relativistic plasma in the 
radio lobes (77) would have to be fc/77 ~ 1.3 x lO'^ if 
no other form of pressure support is present. Simi- 
lar (although somewhat smaller) discrepancies in pres- 
sures are seen in several cooling core sys tems such as 
Hydr a A ijMcNamara et al.ll2000r. Perseus llFabian et all 
12003), and Abell 2052 l|Blanton et al.ll200i| . These ob- 
servations suggest that the presence of an additional 
form of pressure support (such as a hot, thermal com- 
ponent) within the cavities may be necessary for these 
systems. The isobaric radiative cooling time of the fil- 
ament is Tcooi ~ 8.4 X 10^ yr. This cooling time is 
significantly shorter than the age of the cluster, but is 
much longer than the synchr otron age of the rad io tail 
(rtaii = 1.1 X 10^ yr) found bv lTavlor et all fl99l . This 
suggests that the gas in the X-ray filament did most of 
its cooling to the current temperature in the dense clus- 
ter core and was later displaced to the current position 
through the influence of the radio source. We determine 
a mass of the filament of ^ 1.5 x 10^° Mq which gives a 
cooling rate in the filament of ~ 2.5 Mq yr~^. 

At high resolution, the inner regions of the radio source 
PKS 1508-1-059 are resolved into a compact core and two 
bright knots located approximately 1'.'3 on either side 
of the core. The narrow inner region appears to de- 
collimate at a distance of 2" (~ 3 kpc) along each of 
the oppositely directed radio jets. The location of the 



de-coUimation is co-incident with a sharp drop in the X- 
ray surface brightness as seen in Figure El The X-ray 
core appears to have an hour-glass shape with the radio 
jets propagating along the pinch axis. The sharp X-ray 
surface brightness drop suggests that the de-coUimation 
of the radio jets is likely due to a decrease in the confining 
pressure of the surrounding medium. I n a study of the 
flaring of the inner radio jets of 3C 31, ILaing fc Bridle! 
()2002t) find that the radio jets are overpressured by a 
factor of 8 at the beginning of the flaring region. 
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Fig. 6.— The central 25" X 25" (36 X 36 kpc) of the adaptively 
smoothed Chandra image shown in Figure The contours show 
the 8350 MHz radio emission of PKS 1508-(-059 from Taylor et a^ 
11994*) . The X-ray emission shows the broad hour-glass shaped 
cluster core. The bright radio core is co-incident with the brightest 
region of the X-ray core. The radio jets propagate along the pinch 
axis of the core and are well coUimated on either side of the core 
out to roughly 2" (~ 3 kpc). Beyond this distance, the radio jets 
flare outward at the sharp drop in X-ray surface brightness. 

Further out, they find that the jets of 3C 31 return to 
roughly pressure equilibrium with the external medium. 
An alternative method to de-coUimate the jets would be 
through the interaction of the radio jet with dense gas 
clumps in the ICM. Unfortunately there is insufficient X- 
ray data in the jet and de-coUimation regions to allow us 
to determine the pressure and spectral properties in the 
surrounding thermal medium. We will investigate this 
connection in more detail with our deep (80 ksec) Cycle 
5 Chandra data set. Beyond the flaring point, the 8.4 
GHz radio jets continue to propagate to a distance of ~ 
10" (~ 15 kpc). The northern jet appears to follow a 
linear trajectory, while the tip of the southern jet shows 
a curvature toward the south. The radio morphology 
of the southern jet appears to be traced by the cool X- 
ray filament which extends south of the cluster core (see 
S I5.2|l . The radio spectral index [Si, cx j/") of the inner 
jets of PKS 1508-1-059 runs from a = -1-0.21 in the core to 
'--^ a = —2.7 at the end s of the jets betwee n frequencies 
of 4860 and 8515 MHz l|Tavlor et al.lll994j) . This steep- 
ening of the spectrum along the jets is consistent with 
the confinement o f the synchrotron pl asma by the exter- 
nal thermal ICM. iTavlor et aP l|1994j) also find that the 
radio spectrum of the jets shows a curvature such that 
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the spectral index is steeper at higher frequencies. Using 
synchrotron aging arguments and the spectral shape of 
the radio emission, they estimate an age of 9.6 x 10^ yr for 
the inner radio jets, and 1.1 x 10'' yr for the outer lobes, 
where we have converted the lifetimes to our cosmology. 

5.1. Is There X-ray Evidence of a Central AGN? 

Examination of the raw X-ray data shows that there 
is no obvious X-ray point source in the core of Abell 
2029, however, the presence of the central radio source 
PKS 1508-1-059 has led us to search for spectral evidence 
of the central active galactic nucleus (AGN). The cluster 
core displays a broad central plateau which has a ra- 
dius of approximately 6". We have extracted the X-ray 
spectrum from a circular region (r=3", 4.3 kpc) centered 
on the flat-spectrum radio core (|Tavlor et al. 1994). We 
have employed two different techniques to account for the 
background: 1) a local background determined from an 
annulus centered on the core with an inner radius of 3" 
and outer radius of 5", and 2) a blank-sky background 
for the circular region determined from the Chandra Cal- 
ibration Database. After background subtraction there 
were 540 counts in the region centered on the cluster core 
using the local background file and 1970 counts using the 
blank-sky background. The data sets were binned to pro- 
vide at least 25 counts/bin, and a variety of models were 
examined for each data set. In both cases, the best fit 
was obtained for a single temperature MEKAL model, 
and in no case were we able to find a good fit to the 
core region which included a reasonable power-law com- 
ponent. All fits including a power-law component gave 
a steep exponent of T > 3.4, while typical radio-loud 
AGN have T ~ 1.5 l|Mushotzkv et alJTl99^ . We note 
also that there is no obvious point source seen in a hard- 
ness ratio map of the broad cluster core. The broad core 
is consistent wit h the radia l surfa ce brightness profile of 
the cluster bv fLewis et ahl ll2002|) . The lack of spectral 
evidence in the Chandra data for a central power-law 
component suggests that the central AGN may be very 
weak and/or very heavily obscured. 

5.2. Central Temperature Structure 

We have created a temperature map of the central re- 
gion of Abell 2029 to investigate the two-dimensional 
temperature structure and interactions with the radio 
source. The map covers the central 95 kpc of Abell 
2029 and has compact sources removed. Based on the 
lack of evidence for excess emission from the central 
AGN (§ 15. 1|) . we have not excluded data from the clus- 
ter core. The temperature map was created within the 
ISIS software environment (Houck & DcNicola 2000) us- 
ing an adaptive binning technique. For each pixel in 
the output map, the spectrum is fit to the smallest sur- 
rounding region which contains at least 1050 counts. The 
smallest extraction region for our temperature map was 
3'.'5 X 3'.'5 in the image center, and the largest region 
was 15'.'7 X 15'.'7 on the image edge. Each spectrum was 
corrected for a background spectrum extracted from the 
the blank-sky background files from the same source re- 
gion. The spectra were binned to contain at least 20 
counts/bin, and were fitted with a MEKAL model. The ab- 
sorption was fixed at Galactic, and the abundance was set 
to 0.45^0. The resulting temperature map is shown in 
Figure [71 This map has been convolved with a Gaussian 



of width 1.5 pixels to smooth out the inter-pixel varia- 
tions. The temperatures in the map run from roughly 3.8 
keV in the core to ~ 9.5 keV, with the average temper- 
ature around 6.8 keV, and typical errors on kT of order 
15 - 20%. 

3.8 5.2 6.5 7.8 9.5 
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Fig. 7. — Temperature map of the central 95 kpc of Abell 2029 
with the color scale shown on top in units of keV. The cluster 
gas is coolest in the center (around 3.8 keV) and shows cool gas 
associated with the filament near the southern radio tail. The 
northern radio tail appears to be embedded in a region of average 
cluster temperature. The temperature structure in the core of the 
cluster is complex, likely due to the superposition of the spiral 
excess feature and the X-ray filaments. 

The temperature map shows a triangular-shaped re- 
gion of cool gas in the cluster core, surrounded by patches 
of cool and hot gas at larger radii. The central jets of 
PKS 1508-1-059 are completely immersed in the cool gas, 
and the X-ray filament tracing the southern radio tail 
also appears to be cool. In fact, there is additional cool 
gas to the south-east of the southern tail, suggesting 
that it may be surrounded by a co ol shell as is found 
in radio-filled X-ray bubbles (e.g., 'Nulsen e t all 120021: 
Fabian ct al. 2000; Blanton ct al. 2003). There is no ev- 
idence in the temperature map of strong X-ray shocks 
associated with the radio lobes, although we cannot rule 
out the transfer of energy from the radio source to the 
thermal gas through weak shocks or pressure waves as 
seen in Perseus (Fabian ct al. 2003). The situation ap- 
pears to be significantly different for the northern ra- 
dio extension where the gas surrounding the radio tail 
is similar to the average cluster temperature. We note, 
however, that the relation of the gas temperature to the 
radio structure in the core of Abell 2029 is complicated 
by the spiral X-ray excess discussed in § El 

6. SPIRAL EXCESS 

In § O we mentioned that the diffuse X-ray emis- 
sion is not symmetrically distributed about the cluster 
core. In order to investigate this further, we compare 
the smoothed X-ray image to a smooth elliptical isopho- 
tal model. We used a Gaussian smoothed (o'=2") im- 
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age (corrected for background and exposure) as input 
to the IRAF/STSDAS task ellipse. The fits allowed 
the ellipticity, position angle, and intensity of elliptical 
isophotes to vary but the ellipse centroids were fixed at 
the compact X-ray core at RA=15 10 56.11, Dec.= +05 
44 41, co-incide nt with the flat-spe ctrum radio core of 
PKS 1508-^059 ijTavlor et al.llT99|) . Compact sources 
were excluded from the fits by setting a clipping thresh- 
old to exclude pixels which deviate more than Scr from 
the annulus mean. The initial values used for the el- 
lipticity and position angle (PA) were 0.2 and 20°, re- 
spectively. The top panel of Figure |H1 shows the best-fit 
ellipticity and position angle for each isophote where the 
radius of the isophotes follows a geometric progression. 
The data were fit from a radius of 3 kpc (~ 2") to the 
maximum outer semi-major radius of 310 kpc (3^6) al- 
lowed by the size of the ACIS S3 chip. The PA is roughly 
constant at 22° from 80 kpc out to 300 kpc, in good 
agreement with t he optical value of 21° determined by 
lUson et all l)199lD . Interior to 80 kpc, the isophotal PA 
increases to ~ 44° at a distance of 14.5 kpc then drops to 
21° at 2.9 kpc, indicating significant structure in the in- 
ner regions of the cluster. The ellipticity of the isophotes 
is fairly constant at 0.26 outside 38 kpc and increases 
to 0.3 around 34 kpc then drops to 0.07 at 2.9 kpc. In 
comparison, the diffuse R-Band e mission of the cD halo 
has an ellipticity of 0.6 l)Uson e t al. 1991). A detailed 
discussion of the intrinsic shape of Abell 2029 based on 
th e observed optical and X-r ay ellipticities is presented 
bvlBuote fc CanizaresI l(T? 
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Fig. 8. — Plots of the position angle and ellipticity of the best fit 
elliptical model to the Gaussian smoothed ((t=2") Chandra image 
of the central region of Abell 2029. Top panel shows fits to the 
full image and bottom panel shows fits after masking the region 
of excess emission. The PA and ellipticity in the bottom image 
are more sensitive to features in the X-ray surface brightness due 
to the smaller averaging areas in each annulus. The spiral galaxy 
absorption feature (Clarke et al. 2004) is located at a distance of 
roughly 130 kpc where the bottom panel shows a variation in the 
PA. 

The fit parameters were used to create a smooth ellip- 
tical model of the emission which was subsequently sub- 



tracted from the Gaussian smoothed X-ray image. The 
resulting residual image (blanked outside the ellipse 
fitted region) is shown in Figure IHl The residuals show 
a striking dipolar (i.e., one-armed) spiral pattern which 
extends from the core in a clockwise direction outward 
to a position angle (east of north) of roughly 225° and 
radius of ~ 1.'5 — 2'. To investigate the possibility of 
absorption causing the feature, we have split the Gaus- 
sian smoothed (cr =2") Chandra X-ray emission into 
soft and hard X-ray bands (0.3 - 1.0, and 2.0 - 10.0 
keV) and created a hardness ratio image of the cen- 
tral region of the cluster. We do not see any evidence 
of hard X-ray emission expected from photoelectric ab- 
sorption near the cluster core. This indicates that the 
feature is likely due to excess X-ray emission. Figure |21 
also clearly reveals the unrelated X-ray absorption 1.'5 
south of the cluster core (seen clearly in our hardness 
ratio image ) due to the disk of a foreground edge-on spi- 
ral galaxy ijClarke et alJl20f3l . The positive and nega- 
tive deviations of the dipolar feature are at the ~ 15% 
level compared to the smooth elliptical model, indicat- 
ing roughly a 30% excess. A morphologically similar spi- 
ral X-ray feature is seen in the Chandra observa tions 
of the Centaurus cluster l|Sanders fc Fabicinll2002() . In 
the case of the X MM- Newton observations of Perseus 
(Churazov. Forman. Jones. & Bohringcr 2003), the ap- 
parent spiral visible in the surface brightness deviation 
image is a result of the interaction of the lobes of the 
central radio source (3C 84) with the dense ICM. This 
does not appear to be the case for Abell 2029 as the ra- 
dio lobes are much narrower than those of 3C 84, and 
the spiral excess extends well beyond the compact radio 
source. 
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Fig. 9. — Residual image of the central region of A2029 after 
subtracting a smooth, elliptical isophote model. The model was 
obtained by fitting the Gaussian smoothed (a=2") image to an 
elliptical model allowing the ellipticity and position angle to be free 
in the task ellipse in IRAF but with the core position held fixed. 
The image displays a dipolar spiral pattern. The linear feature 1.5' 
south of the cluster core is due to p hotoelectric absorpt ion from the 
disk of a foreground spiral galaxy IClarke et aljl2004l^ . 
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We have examined the robustness of the excess emis- 
sion feature by running a variety of fits to the smoothed 
X-ray emission with: 1) the centroid, ellipticity and PA 
free, 2) the centroid and ellipticity fixed and PA free, 3) 
the centroid and PA fixed and ellipticity free. All fits pro- 
duced the same excess emission region with only slight 
variations in the structure in the central region. We have 
also used the model-subtracted image to create a smooth 
masked area blocking out the excess emission feature. 
This mask was then included within the ellipse task 
in IRAF to produce a smooth model of the remainder 
of the cluster emission. This technique has the advan- 
tage of avoiding over-subtraction of the cluster emission 
due to averaging within annuli that include the excess. 
The residual image from the masked fitting technique is 
shown in Figure [TUI and the PA and ellipticity are shown 
in the bottom panel of Figure |H1 
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Fig. 10. — Similar residual imago to Figure |5] but with the 
model obtained by masking the region of excess emission in Figurelol 
before fitting the elliptical model. This technique reveals the same 
excess emission structure but reduces the over-subtraction effects 
seen in the previous image. 

To investigate the details of the dipolar spiral feature 
we have used the positive excess to set an extraction re- 
gion from the core out to a radius of ~ 40" (58 kpc) at a 
position angle of 50°. Using CALDB v2.26, we extracted 
the spectrum of the excess from this region, and also ex- 
tracted a comparison spectrum by rotating the region 
mask by 180° about the cluster center. In the energy 
range of 0.7 — 8.0 keV, the two regions contained 18730, 
and 14610 background-subtracted counts for the positive 
and ambient areas, respectively. The data were fit with a 
single MEKAL model within XSPEC with the absorption 
fixed to Galactic, and the temperature and abundance 
left as f ree parameters. The fit results are shown in Ta- 
ble Ell The region containing the excess is found to be 
cooler (by ^ 0.9 keV) than the ambient region at the 90% 
confidence level. Since much of the emission in these two 
spectra will be due to projected foreground and back- 
ground gas, it is likely that the real local temperature 



contrast is larger than shown in Table IA3I In compar- 
ison, the plume in the Centaurus cluster is found to be 
roughly 0.25 keV cooler tha n the ambient temperature 
at a similar radius ( Sanders fc Fabian 2002). In both the 
iSanders fc Fabianl ( 2002) study of Centaurus cluster and 
our analysis of Abell 2029, the metallicity of the X-ray 
spiral is found to be consistent with that in the ambient 
gas. 

One possibility for the origin of the spiral excess is that 
it is a cooler gas cloud which has been stripped as it fell 
into the center of the cluster. The gas may have been as- 
sociated with a galaxy or group of galaxies which merged 
with Abell 2029, or it might have been purely gaseous, 
or gas occupying a dark matter (DM) potential well. We 
have estimated the mass of gas required to produce the 
observed excess. We assumed that the spiral excess re- 
gion was located in the plane of the sky at the distance 
of the cluster center. We assumed that the thickness of 
the spiral excess region along our line-of-sight was equal 
to the width of the feature in the image (i.e., that it 
is a spiral tube of excess density). The emissivity was 
determined from spectral fits to the region. We depro- 
jected the foreground and background emission to obtain 
the local density in the feature, and integrated this over 
the volume of the region. The estimated total mass is 
-^^spiiai ~ 6 X 10^^ Mq. It is likely that this is a lower 
limit to the mass required. A quick examination of the 
XMM/Newton archive data on Abell 2029 confirms the 
existence of this spiral feature, but indicates that it ex- 
tends about a factor of 2 further out in radius beyond 
what is seen in the Chandra image. The mass in the 
spiral is increasing rapidly with radius over the range 
observed with Chandra, so this suggests the total mass 
might be more than a factor of two larger. The mass 
would also increase if the feature were not located in the 
plane of the sky. With this gaseous mass, the material 
could not have come from a single galaxy; it would have 
to come from a group of galaxies or gas trapped in a dark 
matter potential well. The mass of the X-ray "plume" 
seen in the Centaurus cluster was estimated to be more 
than 1000 times smaller (Sanders & Fabian 2002) and is 
more likely due to a single infalling galactic system. 

In Appendix ^ we show that a spiral trajectory is 
a solution for an object with angular momentum sink- 
ing into the center of a cooling flow cluster. This so- 
lution applies if the object is sufficiently heavy, has a 
nonzero angular momentum, and is falling in at the "ter- 
minal velocity" (i.e., any initial transients have died away 
and drag forces nearly balance gravity) . In polar coordi- 
nates (radius r, angle 9), th e or bit is a logarithmic spiral, 
r = 7-0 exp[a{6 — 9q)] (eqn. IA8(I . A crude fit to the cen- 
ter of the spiral excess in Abell 2029 suggests Bq = 0, 
ro ^ 40", and a ~ —0.4, if 9 is the position angle mea- 
sured from north to east. The constant a determine s th e 
pitch angle of the spiral, and is given by equation l|A7|l . 
Based on gas density from a deprojection of X-ray surface 
brightness in the non-spiral-excess region, the required 
mass column density of the sinking object is E ~ O.OIC 
g cm~^, where C ^ 1 is the drag coefficient of the body. 
This implies that the total mass of the sinking object is 
A/tot ~ 4 X 10^ C (D/10kpc)2 Mq, where D is the diam- 
eter of the object. For D ^ 10 kpc, this mass is much 
smaller than the gaseous mass needed to produce the spi- 
ral excess. The two masses would be similar if D > 100 
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kpc, but this is larger than the inner parts of the spiral 
excess. This suggests that the sinking object was purely 
gaseous, and was being ablated to produce the spiral as 
it fell into the cluster center. Unfortunately, this compli- 
cates the calculation of the dynamics; clearly, numerical 
hydrodynamical simulations would be useful to evaluate 
this model. 

Alternatively, the spiral feature might be the result of 
sloshing motions in the core of the cluster, perhaps in- 
duced by a past cluster merger (Maxim Markevitch, pri- 
vate communication). Sloshing features, which are simi- 
lar to "cold fronts " , have been seen in a nu mber of cooling 
flow clusters fe.g.. lMarkevitch et a l.*20'o7l. but not in the 
form of extended spirals. A spiral sloshing feature might 
arise in an offset merger with significant angular mo- 
mentum. Short spiral features are seen in some merger 
simulations ('e.g.. lRicker fc SaraziiJl200'lL Figure 7), but 
they generally are stronger shocks and occur during very 
strong mergers. Whether weaker spiral sloshing features 
can persist long after a merger in an apparently relaxed 
cluster like Abell 2029 is uncertain; again, numerical hy- 
drodynamical simulations would be useful. 

The central radio source PKS 1508-1-059 reveals a clas- 
sical C-shaped WAT morphology. As discussed in § [S] 
the northern radio extension to PKS 1508-1-059 follows 
roughly the same position angle as the inner jets, but 
is offset ~ 6" (8.6 kpc) west of the main jet axis. The 
outer southern tail is also offset to the west and appears 
to be rotated by roughly 90° from the inner jet axis. 
In standard WAT radio sources, this C-shaped morphol- 
ogy is attributed to relative motion between the radio 
galaxy and the surrounding intergalactic medium. It has 
been argued that the orbital motions of a cluster-core 
cD radio galaxy a r e too small to produce th e distortions 
ijEilek et al.lll984t lO'Donoehue et alJ 11993(1 . In classi- 
cal swept-back cluster-core sources, the morphology is 
therefore often attributed to a merger where the ICM 
is moving past the radio source l|Burns eFaDIHll). In 
Abell 2029, there is no major cluster merger apparent, 
although we do see the spiral excess which is likely due to 
either an infalling system or sloshing of the cluster core; 
thus it is possible that the WAT morphology is due to 
an interaction between the X-ray gas in the spiral fea- 
ture and the central radio source. Figure ITTI shows the 
central 95 kpc of the residual image with the 1490 MHz 
radio contours of Tavlor et al. ( 1994) overlaid. In projec- 
tion, the excess appears to run directly between the inner 
southern jet of PKS 1508-1-059 and the steep-spectrum 
radio tail. Due to the uncertainties in the origin of the 
spiral excess, it is difhcult to determine if the offset be- 
tween the inner and outer southern radio structures is 
due to an interaction of the two features, or if it is sim- 
ply an unrelated projection effect. The radio tail itself 
appears to lie in a local depression, with a very com- 
pact, deep X-ray depression co-incident with the south- 
ern tip of the radio tail. Comparing the 0.3 — 10 keV 
counts in a circle of radius 2" centered on the depression 
to a surrounding annulus between 2" and 12" we find 
that the depression is significant at the 4cr level. Unfor- 
tunately, the small number of counts in the depression 
do not allow us to undertake a more detailed analysis. 
The northern radio lobe appears to be unaffected by the 
spiral feature. If the outer radio lobes represent (undis- 
turbed) buoyantly rising bubbles of relativistic particles 



and magnetic fields, we would expect them to rise along 
the lowest thermal pressure gradient (i.e. along the minor 
axis of the X-ray emission). In this scenario, the most 
likely source configuration would be an S-shaped source 
with the northern lobe rising to the north-west and the 
southern lobe rising to the south-east. 
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Fig. 11. — Residual image of the central 95 kpc of Figurc lT?)l with 
VLA 1490 MHz radio data of ITavlor eFall 119941) overlaid. The 
excess emission region runs (in projection) between the southern 
jet and outer radio lobe. 

7. SUMMARY AND DISCUSSION 

We have presented a detailed analysis of the Chandra 
ACIS-S3 observations of the central regions Abell 2029. 
Although this cluster appears to be remarkably relaxed 
on large scales, there is significant structure within the 
central 200 kpc radius. In general, the X-ray emission 
is elliptical along a position angle of 22°, consistent with 
the position angle of the optical emission of the central 
cD galaxy. The thermal emission is asymmetric about 
the cluster core with excess emission to the north-east 
and south-east compared to the south-west and north- 
west, respectively. The extent of the excess emission is 
seen more clearly after we have fitted and subtracted a 
smooth elliptical model from the data. The resulting 
X-ray residuals show a clear one-armed spiral pattern 
running clockwise from the core out to radii of ~ 170 
kpc. The excess X-ray emission is ^ 30% higher than 
the smooth elliptical model. Based on spectral fits to the 
region, we estimate a total mass of the spiral excess of 
Mspirai ~ 6x 10^^ Mq, although this should be considered 
a lower limit due to uncertainties in the size and location 
of the region. It seems likely that the excess feature may 
be the result of stripping of gas from a galaxy group or 
gas trapped in a dark matter potential. An alternative 
possibility is that the excess is the result of sloshing of 
the cluster core due to a past cluster merger. It would be 
helpful to have numerical simulations of both models to 
help assess the feasibility of creating the spiral structure 
from the different mechanisms. 

The central radio source in Abell 2029 displays a 
wide-angle-tailed morphology which is typical of clusters 
undergoing mergers. Several authors have previously 
remarked on the unusual occurrence of such a source in 
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the core of a relaxed clu s ter su c h as Abell 2029 (see e.g. 
T.ewis. Stocke. Ruotel 1200^ IHaTdr.a,st1e Sa,ke1lio.i 
2004 iLecavelier des Etanes. Gopal-Krishna. fc Durret 
20041) . One possible explanation for the WAT morphol- 
ogy is that the mechanism responsible for producing 
the X-ray spiral excess has also influenced the radio 
source structure. The tip of the southern radio jet of 
PKS 1508-1-059 bends toward the south at roughly the 
location where it intersects (in projection) with the inner 
regions of the X-ray excess. In addition, the southern 
radio lobe is displaced (again in projection) to the west 
of the minor axis of the X-ray emission. If the emission is 
due to a buoyantly rising detached radio lobe we would 
expect it to follow the steepest X-ray pressure gradient 
which should be along the minor axis. In comparison, 
the northern outer lobe, which is apparently unaffected 
by the spiral excess, appears to be roughly directed 
along the minor axis. The southern radio lobe is traced 
on the western side by a bright X-ray filament which 
may be a partial shell similar to that seen for Hydra A 
|McNamara ct al. 2000) . The temperature map of the 
central cluster region shows that the gas surrounding 
the southern radio lobe is cool, consistent with the 
results found for the brig ht rims around radio bubbles in 
clusters such as Perseus (*F abian et a l.' 2000") and Abell 
2052 (Blanton et al. 2001). We will combine the 19 ksec 
archival data with our deeper 80 ksec Cycle 5 data set 
to examine the region of the radio foot in more detail to 
determine if there is a complete rim. 

We find that the total cluster spectrum of Abell 2029 
is poorly fit by a single temperature thermal model. The 
best fit is obtained for a two temperature MEKAL model 
with high and low temperatures of kT^igh — lAT^'"^ 
keV and kTiow — 0.111;q3 keV, and an abundance of 
Z = 0.49t;o6 Zq. a nearly identical fit quality is ob- 
tained for a cooling flow model with a similar metal- 
licity, gas cooling over a range of temperatures from 
7.94 keV down to the hmits of the model at 0.08 keV, 
and a mass deposition rate of M — 56^.21 yr~^. 
The range of temperatures from the cooling flow model 
is more than a factor of 50 between the hot and cool 
components. This large temperature range is relatively 
unique as most cooling flow clusters show gas cooling 
down to only ~ 1/2 — 1/3 of the ambient temperature 
(jPeterson et al.ll2'003(l . A recent study of the central 20 
kpc radius of Abell 2029 with FUSE found no evidence 
of the O VI doublet emission which is a good tracer o f 
gas at 10^ - 10^ K l|Lecavelier des Etangs et all 12004"). 
Based on this non-detection, they place upper limits on 
the cooling flow rate through 10^ K of 13 Mq yr~^. This 
limit is lower than the 56 Mq yr~^ we determine to be 
cooling through ~ 10^ K based on the Chandra data. 
The discrepancy in cooling rates may be largely due to 
the fact that our cooling rate is determined over a much 
larger aperture of radius 167 kpc. In a future paper, we 
will use the deep combined Chandra data to investigate 
the spectral properties in more detail. 

Abell 2029 may be a very young cooling core cluster 
where the gas has only recently started (or re-started) 
cooling to low temperatures. In this scenario, the radio 
source may have only recently re-started as a result of 
feeding from the cooling flow, and thus feedback mech- 
anisms may not have had time to suppress the cooling 



to low temperatures. The cooling flow may have only 
existed for a short time, and too little gas has cooled 
to low temperatures to produce optical emission lines or 
star formation. One concern with this feedback scenario 
is that gas must have cooled and reached the central 
black hole in order to produce the current radio core. It 
may be that gas in Abell 2029 is presently cooling for 
the first time; alternatively, a previous cooling episode 
may have been disrupted by a merger or heated by the 
central AGN activity. The presence of the spiral excess 
feature may suggest that Abell 2029 has had a minor 
merger in the not too d i stant past, and numerical simu- 
lations bv lGomez et alJ l|2002t ) indicate that mergers can 
lead to the disruption of cooling flows in clusters. One 
argument against a recent merger or hydrodynamical dis- 
ruption of the core of Abell 2029 is the stron g increase in 
the Fe a bundance in the c l uster core seen bv lLewis et alJ 
()2002|) . iBohringer et al.l l)2004j) argue that mergers or 
other processes which disrupt the cooling core might also 
mix the gas and hence eliminate such abundance gradi- 
ents. In the case of the minor merger in Abell 2029, the 
core may not have been signiflcantly disrupted and/or 
the infalling material may have had an enriched metal- 
licity which could contribute to the observed enhance- 
ment. Unfortunately spectral flts to the current data do 
not provide useful limits on the metallicity of the excess 
material relative to that of the ambient gas. In addition, 
it is possible that the onset of cooling was delayed by 
an AGN heating process which did not mix the cluster- 
center gas. 

Our analysis of the Chandra data reveal the presence 
of significant X-ray structure as well as multiphase gas in 
the core of Abell 2029. This result is in agreement with 
previous X-ray studies of the cluster based on ROSAT 
and ASCA data fSaraz in et ail 119921: lEdee et all [1991 
iPeres et al. 1998; Saraz in et al.lll998D but is inconsistent 
with the c onclusions of t he Chandra analysis of the clus- 
ter bv .Lewis et al.l l|2002|) who found no evide nce of multi- 
phase gas or significant X-ray substructure. iLewis et alJ 
()2002t) fit the cooling flow model individually to each of 
the inner three annuli (roughly inner 35 kpc). Unfortu- 
nately, at the time of their analysis, the QE degradation 
of Chandra at soft energies was not known, thus their 
spectral models would over-predict the flux at low ener- 
gies resulting in weaker constraints on the cool compo- 
nent of the emission. The presence of some X-ray struc- 
ture in the central region of Abell 2029 was also noted by 
ILewis et a1 .1 lf2?i?il : in a Gaussian (ct = 2") smoothed im- 
age they saw position-angle-dependent variations in the 
surface brightness within the central 1' region where we 
find t he spiral excess feature. Lewis, Buotc. & Stock3 
(|2003) used the Chandra data to measure the density 
and temperature distribution in Abell 2029 down to the 
inner 3". Based on the apparent lack of both multiphase 
gas and signiflcant X-ray structure they used the assump- 
tion of hydrostatic equilibrium to determine the distri- 
bution of total mass with radius in the system. Fitting 
the mass profile to various m odels they found that the 
distribution was well fit by the lNavarro. Frenk. fc Whitd 
llT99l (NFW) profile over the range of radii between 5 
and 260 kpc. Th e lack of deviation s of their data from the 
NFW profile led ILewis et all ('2003) to conclude that the 
fiatte ned dark matter profile pred icted for self-interacting 
DM ijSpergel fc Steinhardti I2fl0(1) was inconsistent with 
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the Abell 2029 data. The presence of significant X-ray 
structure and multiphase gas that we find within the 
central 130 kpc indicates that the central regions of the 
cluster are not fully in hydrostatic equilibrium. Without 
detailed numerical simulations of this system, it is diffi- 
cult to say at what level the d eparture frorn hydr ostatic 
equilibrium has influenced the iLewis et all l|2003(l fits to 
the mass profile. 
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APPENDIX 

SPIRAL SOLUTION FOR SINKING BODY 

Here, we derive the equations for the orbit of a sinking body in a cluster of galaxies subject to gas drag and gravity 
forces. We also present an illustrative solution which has the form of a spiral. We will assume that the cluster is 
spherical and static. Then, by symmetry, the orbit will be planar. We use the polar coordinates r and 6 to describe 
the orbit, where r is the radius out from the cluster center, and 9 is the angle around the orbital plane, measured from 
any arbitrary axis. Let p{r) be the ambient gas density in the cluster. Assume that the mass m and cross-sectional 
area A of the sinking body are constant. Let v be the velocity vector of the body. We will write the drag force on 
the sinking body as Fdrag ~ ^(1/2) p{r) AC{v'^) wv, where C(w^) is the drag coefficient C. The speed of the body is 



given by u = Y -|- {rO)"^ , where the dots denote time derivatives. Then, the equation for the angular momentum of 
the sinking body is 

d{r^e) I p{r)C[P + {ref] 



dt 2 Y. 



where E = m/A is the mass column density of the body. Similarly, the equation for the radial acceleration is 

.. -2 GMjr) 1 pir)C[r^ + ir9n / ^ , . 

r - rO = ^ \ r"^ + (rO)'' r , (A2) 

r^ 2 2j ^ 

where M{r) is the total clus ter mass interior to r. We assume that m ^ M(r) at all radii of interest. The second and 
third terms in equation (jA2|l are the centrifugal and gravitational accelerations. The final terms in both equations (jAl|l 
& ljA2|l represent the effects of drag. 

For our simple illustrative solution, we will assume that the drag coefficient C is constant. We will adopt a singular 
isothermal sphere model for the mass distribution in the cluster, which implies that M{r) = 2a'^r/G, where a is the 
one-dimensional velocity dispersion in the cluster. Since we are mainly interested in the central regions of Abell 2029, 
which is a cooling core cluster, we adapt a simple form for the intracluster gas density distribution, p{r) = p{){r^/r). 
This distribution is roughly consisten t with the observed gas distributions in the centers of many cooling flow clusters 
(e.g.. IWhite. Jones, fc FormanllT99'^ . Here, we will use "0" subscripts to denote the initial position of the sinking 
body. With these assumptions, the radial and angular equations become 

2^2 1 PoroC , / ■ r 



rO' ^-—--^^^^f^ + irOf-, (A3) 



and 

djr-'e) _ 1 porpC 
dt 2 1] 



+ {reYre. (A4) 



Any detailed discussion of the full set of solutions to equations ljA3|l & ljA4p is beyond the scope of this observational 
paper. However, we note that one set of solutions is given by 



and 



r = ro (l + Ooat) , (A5) 
= i In (l + e^at) + 00 . (A6) 
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Here, tq and Oq are the initial radius and azimuthal angle at t — 0, and Oq is the initial value for the angular velocity. 
The constant a is given by 

_ sign(^o) 



\poroC J 



(A7) 



- 1 



The geometric shape of the orbit in this solution is that of a logarithmic spiral, 

r = ro exp[a(6' - 6*0)] . 



(A8) 

Equations IjASp . (|A6|) . & IjASp do not represent a general solution of the equations. They only apply exactly when 
the initial radial and angular velocity of the sinking object are given by 



ro = — 



c rporpC 



and 



= ± 



\ 



2cr2 



V 2S 



(A9) 



(AlO) 



However, one can show that this is an asymptotic ("terminal velocity") solution for a sinking body with angular 
momentum which applies after transients associated with the initial conditions have died away. The solution also 
requires that the mass column of the sinking body exceed a lower limit, S > po^oC'/2. When this limit is approached, 
the motion becomes radial infall {9 = 9q — constant) at the terminal velocity for radial motion r = ro — \/2(7t. For 
radial motion, the terminal velocity is \/2a under the assumptions made here. 
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Table A1. XSPEC Fits to the Inner 116" Radius Region 



Model 


Nh 




fcThigh 


Abundance 


M 


xVd.o.l. 




(xlO^o cm-2) 


(keV) 


(kcV) 


(Solar) 


[Mq yr-i) 





MEKAL 


(3.14) 




7 Q7+0. IG 
-0.14 


J-' -0.06 




532/432= 


1 


23 


MEKAL 


1 97+0.64 
-0.43 




7 7(1+0.35 
'•'^-0.13 


49+0-03 




508/431= 


1 


18 


MEKAL.e 


(3.14) 




'•^*-0.12 


47+0.03 
^•^'-0.05 




499/411= 


1 


21 


MEKAL.e 




'^■^^-0.03 


'■^"-0.17 


p;n+0.04 




473/410= 


1 


15 


MEKAL+MEKAL 


(3.14) 


'•^'-0.15 


49+0.03 
'^•^^-0.06 




507/430= 


1 


18 


MEKAL+MEKAL 


— .75 




7 fil+o-^" 


50+0-02 




506/429= 


1 


18 


MEKAL+MEKAL.c^ 


(3.14) 


O.14l0;0i 


7 47+0.23 
-0.14 


n 50+0-03 

*J-''^-0.07 




472/409= 


1 


16 


MEKAL+MKCFLOW 


(3.14) 






50+0-03 


561^1 


513/430= 


1 


19 


MEKAL+MKCFLOW 




(0.08) 


7 oQ+0.26 


U.01_g Q5 




509/430= 


1 


18 


MEKAL+MKCFLOW.C 


(3.14) 




o ,-,0+0.50 


51+0-0-1 
'-'-^1-0.06 


62+1® 
°^-21 


477/409= 


1 


17 


MEKAL+MKCFLOW.e 


n QQ+-''8 

u.»»_Q 52 


(0.08) 


7.78l°;?g 


51+0-01 


9+9 


472/409= 


1 


15 



Note. — Values in parentheses were held fixed in the models. Models -with js are fit over the 0.7 — 8.0 keV range excluding the 1.8 — 2.1 keV 
energy interval, while the remainder of the models cover the entire 0.7 — 8.0 keV range. 

''A two temperature MEKAL model with free absorption excluding the 1.8 — 2.1 keV region resulted in an unconstrained temperature on 
the second component and is not shown here. 



Table A2. XSPEC Fits to the Outer 116" < r < 174" Region 



Model 


Nh 




^^high 


Abundance 


M 


xVd.o.f 






(Xl020 cm-2) 


(keV) 


(keV) 


(Solar) 


{Mq yr-1) 






MEKAL'' 


(3.14) 




~ Q7+0.54 
'•="-0.36 


„ 04+0.O8 




436/343=1 


27 


MEKAL.e^ 


(3.14) 




7 00+0-54 
'•y^-0.34 






406/322=1 


26 


MEKAL+MEKAL 


(3.14) 


n+0-04 
'J-1'^-0.03 


o 71+0.57 
°-'l-0.41 


42+0-09 
^•^^-0.11 




399/341=1 


17 


MEKAL+MEKAL 




14+0-02 
"■1^-0.01 


o 71+0.58 
°-' 1-0.39 


40+0-11 
'-'•^^-0.09 




399/340=1 


17 


MEKAL+MEKAL.e 


(3.14) 


n 14+0.03 

'^•l+-0.03 


o 71+0.57 
°-'l-0.39 






363/320=1 


13 


MEKAL+MEKAL.C 


0.84lJ-ji 


n 14+0-04 

'J-1+-0.05 


9.29lO;It 


0.4410.10 




362/319=1 


14 


MEKAL+MKCFLOW'' 


(3.14) 


n 09+0-30 
'J-'jy-o.oi 


10.8711^6 


0.40l°;l? 


27+1 


406/341=1 


19 


MEKAL+MKCFLOW.e 


(3.14) 


0.081^0^ 


ll.8l±l;- 


o.4iiri§ 


331^ 


366/320=1 


14 


MEKAL+MKCFLOW.e 






10.5911:^1 




12+f 


362/319=1 


13 



Note. — Values in parentheses were held fixed in the models. Models with _c are fit over the 0.7 — 8.0 keV range excluding the 1.8 — 2.1 keV 
energy interval, while the remainder of the models cover the entire 0.7 — 8.0 keV range. 

''The single temperature MEKAL models with free absorption resulted in absorbing columns consistent with zero for both energy range fits 
and thus they are not shown here. 

''The MEKAL + MKCFLOW model for the 0.7 — 8.0 keV energy range with free absorption was not able to constrain kTi^^ and is not shown 
here. 



Table A3. Spectral Fits to the Spiral Excess Feature 



Region 


kT 
(keV) 


Abundance 
(Solar) 


xVd-o.f. 


Excess 
Ambient 


5.53 23 

6.4i+;fo 


0.7ll;ii 


272/229=1.19 
234/207=1.12 



